The high energy spectrum of NGC 4151 
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ABSTRACT 

We present first INTEGRAL observations of the type 1.5 Seyfert galaxy NGC 4151. Com- 
bining several INTEGRAL observations performed during 2003, totaling ~ 400ksec of exposure 
time, allow us to study the spectrum in the 3 - 300 keV range. The measurements presented 
here reveal an overall spectrum from X-rays up to the soft gamma-rays that can be described by 
an absorbed (Ah - 5 x 10 22 cm -2 ) and non- variable thermal component, plus a Fe Ka line, and 
an exponential cutoff occurs at 110 keV, consistent with earlier claims. The Galactic hydrogen 
column density in the line of sight is Ah .Gal — 2.1 x 10 2o cm“ 2 . The time resolved analysis shows 
little variation of the spectral parameters. The comparison with CGRO/ OSSE data shows that 
the same spectral model can be applied over a time span of 15 years, while the flux varied by a 
factor of 2. Applying a Compton reflection component improves the model fit to the INTEGRAL 
data. Nonetheless the data available to date cannot significantly confirm or exclude the existence 
of reflection, nor is a high iron overabundance in the absorber, as had been previously suggested, 
clearly detectable. 

Subject headings: galaxies: active — galaxies: individual (NGC 4151) — gamma rays: observations — 
X-rays: galaxies — galaxies: Seyfert 
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1. Introduction 

The Seyfert galaxy NGC 4151 is one of the 
most intensly studied AGN (Active Galactic Nu- 
cleus). This is due to the fact that it is one of 
the apparently brightest (V ~ 11 mag, bolometric 
luminosity ~ 10 44 ergs _1 ) and nearest (16.5 Mpc 
for Ho = 75kms“ 1 ) among its class. Observations 
throughout the electromagnetic spectrum revealed 
a complex and somewhat special case of a Seyfert 
1.5 galaxy (see e.g. Dermer & Gehrels (1995) and 
Ulrich (2000) for a general review). In the optical 
the host galaxy appears to be a large barred spi- 
ral with nearly face-on orientation, showing two 
thin spiral arms (Arp 1977). In the high-energy 
domain, NGC 4151 has been studied by various 
X-ray missions. Here the AGN shows a relatively 
hard spectrum with a photon index T which varies 
between 1.4 and 1.7. The spectrum seems to be 
best described by a power law with an exponen- 
tial cut-off at > 50keV and an iron fluorescence 
Ka line at 6.4 keV (Piro 1999; Piro et al. 2000). 
In addition to the fluorescence line, which re- 
veals the existence of cooler material surround- 
ing the AGN, a reflection hump in the 10 - 
30keV has been reported in Ginga/OSSE data 
(Zdziarski et al. 1996) and in BeppoSAX data 
(Petrucci et al. 2001; Schurch & Warwick 2002). 

In this paper we present analysis of recent ob- 
servations of NGC 4151 performed by the INTE- 
GRAL satellite, and compare the results with pre- 
vious studies. 

2. Simultaneous optical, X-ray and gamma- 
ray observations 

Observations in the X-ray to soft gamma-ray 
domain have been performed in late May 2003 by 
the INTEGRAL satellite (Winkler et al. 2003). 
This mission offers the unique possibility to per- 
form simultaneous observations over the 2 - 
8000 keV energy region. This is achieved by the X- 
ray monitor (2-30 keV) JEM-X (Lund et al. 2003), 
the soft gamma-ray imager (20-1000 keV) ISGRI 
(Lebrun et al. 2003), and the spectrograph SPI 
(Vedrenne et al. 2003), which operates in the 20 
- 8000 keV region. In addition to these data an op- 
tical monitor (OMC, Mas-Hesse et al. 2003) pro- 
vides photometric measurements in the V band. 

In order to extend the spectral coverage, we use 
data from RXTE. 
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2.1. INTEGRAL data 

For the analysis of the high energy spectrum 
of NGC 4151 we used the data which were taken 
within 10 degrees radius around the source. The 
journal of observations is shown in Tab. 1. The ex- 
posure times were taken from the ISGRI effective 
on-source times. This value is approximately the 
same for the spectrograph SPI, but for the JEM-X 
and OMC monitors it has to be taken into account 
that they cover a much smaller sky area. Thus in 
the case of dithering observation, the source is not 
always in the field of view of the monitors. The 
analysis was performed using version 4.2 of the Of- 
fline Science Analysis (OSA) software distributed 
by the ISDC (Courvoisier et al. 2003b). 

NGC 4151 was too faint above ~ 150 keV to 
be seen by IBIS/PICsIT. Summing all available 
INTEGRAL data together shows that NGC 4151 
is detected within the exposure time of 408 ksec 
up to 300 keV. No intercalibration had to ap- 
plied between JEM-X and ISGRI, but the SPI 
data appear to be at a higher level by a factor 
of 1.3. The SPI data further show some energy 
bins where the reconstruction software SPIROS 
(Skinner & Connell 2003) could not find a solu- 
tion, and others which indicate problems in the 
background subtraction. This is not surprising as 
most of the data have been taken in staring mode, 
which is known to cause problems in the SPI anal- 
ysis, because the background level cannot be de- 
termined (Dubath et al. 2005). 

The ISGRI spectra have been extracted using 
the standard method. Concerning the calibration 
of the data it has to be taken into account that 
the response of ISGRI is still not well known, even 
though this would be more severe for off-axis ob- 
servations. The spectra of the Crab taken with 
JEM-X show systematic features in the 5 - 7keV 
energy range, which might affect the measurement 
of the iron Ka line. 

2.2. RXTE data 

The All Sky Monitor (ASM) on board the Rossi 
X-ray Timing Explorer ( RXTE) scans about 80% 
of the sky every orbit. This offers a unique way 
to monitor the emission of bright X-ray sources 
like NGC 4151 in the 1.5 — 12keV energy range. 
We extracted fluxes from the RXTE / ASM data 
base. The fluxes have been averaged over one day 


(see the triangles in Fig. 1). The data show a 10 dh 
3 mCrab flux over the duration of the INTEGRAL 
observations. 

2.3. CGRO data 

The source was also monitored in the ^ 
20 — 100 keV band by the CGRO / BATSE ex- 
periment from 1991-2000 (Parsons et al. 1998; 
Harmon et al. 2004). Those measurements are 
made using an Earth-occultation technique, and 
are thus sampled one time per 90 min CGRO or- 
bit. In practice, NGC 4151 is a faint source for 
that instrument, but daily summations lead to 
~ 3 cr determinations. We obtained the BATSE 
data spanning the ~ 3.5 yr overlap of the RXTE 
and CGRO missions, for the purpose of monitor- 
ing the long-term hard-to-soft spectral behavior. 

3. X-ray to Gamma-ray spectrum 
3.1. Time Averaged Spectral Analysis 

The overall spectrum extracted from the INTE- 
GRAL JEM-X and ISGRI data is shown in Fig. 2. 
A systematic error of 3% has been added to the 
ISGRI and SPI data. The fit results however do 
not depend on the SPI data, which do not influ- 
ence any model fit because of the larger error bars 
compared to the ISGRI data in the same energy 
region. ' 

The spectral fitting was done using version 
11.3.1 of XSPEC (Arnaud 1996). The data in 
the range 3 — 300 keV are well represented by 
a rather simple model. We applied a fit of an 
absorbed power law with an exponential high- 
energy cut-off, and added a Gaussian line in or- 
der to account for the iron Ka fluorescence line. 
The best fit values are F — 1.52 ^ ^, Ec = 
109^2 keV> and A H = 5.08+° °® x 10 22 cm~ 2 . 
The fluorescence line shows an equivalent width 
of EW = 119 eV and a flux of f Ka = 4.5 j x 
10 -4 ph cm -2 s _1 (all errors are 90% confidence 
values). The statistics do not allow us to see any 
effects of relativistic broadening of the iron line 
(see e.g. Yaqoob et al. 1995). The source had an 
average flux of /20-200keV = 0.012 ph cm" 2 s -1 = 
9.6 x 10“ lo ergcm“ 2 s“ 1 . This model gives xl = 
1.46 for 189 degrees of freedom. The model of an 
absorbed power law without a cut-off leads to a 
significantly poorer fit result (xl — 1-9). 


Using a more complicate model by adding re- 
flection from cold material to the cut-off power- 
law (the so-called PEXRAV model; Magdziarz & 
Zdziarski 1995) improves the fit. From radio mea- 
surements the inclination angle of the accretion 
disk is estimated to be i ~ 65° (Pedlar et al. 1993). 
The free parameters in the PEXRAV model 
then are best fitted with a photon index of 
T = 1.86 ^o‘o9> a hydrogen column density of 
Ah = 7.0+0® x 10 22 cm -2 , a folding energy of 
Ef 0 id = 414+17® keV and a relative reflection 
of R = 1.04+q 3i . The large uncertainty and 
high value of the folding energy outside the en- 
ergy range covered by the spectrum indicates that 
the reflection component in total is not well con- 
strained. The PEXRAV model leads to xl = 1*19 
for 188 degrees of freedom. 

Another physically meaningful model describ- 
ing Comptonization of soft photons by a hot 
plasma, the so called compTT model, has been 
developed by Titarchuk (1994). The model in- 
cludes the plasma temperature T e of the hot 
corona, the optical depth r of this plasma, and 
the temperature To of the soft photon spec- 
trum. Because the spectrum starts at 2 keV, 
To is not well constrained by the data and has 
been fixed to 10 eV. The model fit resulted in 
Ah = 6.5+q 7 x 10 22 cm -2 , kT e = 35.8+®;^ keV, 
and r = 1.41+0 ^ for plane geometry. The 
compTT model gives xl = 1-38 for 189 degrees 
of freedom, which is significantly worse than the 
PEXRAV model. 

3.2. Spectral Evolution 

In order to study the evolution of the spectrum 
in time, we applied ten time bins over the total ob- 
servation period, combining about 33ksec of effec- 
tive on-time in each bin. We used only those data 
for which we have simultaneous JEM-X and IS- 
GRI coverage. We then applied the cutoff power- 
law model to the single time bins and extracted 
fluxes in the 5 — 20keV and 20 — 100 keV energy 
bin. The results are shown in Fig. 1. No signifi- 
cant variation of the flux can be seen here, and also 
the fit parameters are consistent over the observed 
period within the error bars. In order to look for 
variability on shorter time scales we then binned 
the INTEGRAL JEM-X2 and ISGRI data into 
2000 sec time intervals (Fig. 3). Both lightcurves 
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show a minor outburst at t ^ 1243.0 IJD 1 . An 
analysis of the enclosing time bin 1242.5 IJD < 
t < 1243.5 IJD does show a slightly increased 
flux (/ 20— 200 keV = 0.0132 ph cm” 2 s" 1 = 1.03 x 
10 -9 erg cm -2 s _1 ) but the same spectral parame- 
ters. 

For the same time intervals as used in the high- 
energy analysis, optical photometry values from 
the OMC have been extracted. Due to the high 
signal-to-noise these data have been binned into 
1000 s time intervals (Fig. 4). The minimum in 
the optical lightcurve around t ~ 1243.0 coincides 
with the maximum in the X-ray lightcurve. Be- 
cause of the short time coverage of the observation 
it is not clear whether this is a coniricidence or a 
physical correlation between the optical and the 
hard X-ray domain. 

3.3. Search for 7 -ray line emission 

The SPI spectrometer offers the opportunity 
to look for narrow line features. In case the 
high energy spectrum would be dominated by 
non-thermal processes, one could expect a sig- 
nificant amount of e ± pair annihilation processes 
around 511 keV. At this energy, SPI has a reso- 
lution of FWHM = 1.95 keV (Attie et ah 2003). 
CGROj OSSE was able to determine the 3a up- 
per limit of a broad (450 - 600 keV) emission line 
feature with / < 6 x 10 -5 ph cm" 2 s -1 . We there- 
fore looked for a comparably narrow feature by 
extracting an image in the 499 - 519 keV en- 
ergy band, centred on the 511 keV rest-frame 
energy at NGC 4151’s redshift. No emission 
at the restframe energy of the annihilation line 
is detectable. The 3cr upper limit extracted 
from the SPI data by SPIROS is / 499 - 519 keV = 
10 -4 ph cm -2 s“ 1 . The 3 a line sensitivity derived 
from ground calibration gives a similar value of 
6.5 x 10“ 5 phcm~ 2 s -1 (Roques et ah 2003). The 
measurement presented here is therefore not im- 
proving the CGRO/ OSSE upper limit, because of 
the much shorter INTEGRAL/ SPI exposure time. 

4. Discussion 

The INTEGRAL overall spectrum can be mod- 
elled by a rather simple absorbed cut-off power law 

1 IJ D is related to the modified Juian Date by IJD = 
AMD - 51544.0 


plus an iron Ka fluorescence line. The more phys- 
ical meaningful model of Comptonization of soft 
photons in a hot plasma gives a slightly better fit 
statistic, and shows that those circumstances are 
applicable to the case of NGC 4151. 

The best comparable coverage in the energy 
range studied by INTEGRAL ISGRI and SPI has 
been provided by CGROj OSSE (Johnson et al. 1997). 
The measurements from INTEGRAL and OSSE 
are compared in Tab. 2. The OSSE results have 
been derived by fitting simultaneously the individ- 
ual OSSE spectra, allowing for different normalisa- 
tion. Our results are similar to what has been re- 
ported in Johnson et al. (1997). We compared the 
OSSE data directly with the INTEGRAL observa- 
tion through simultaneous fitting. Figure 5 shows 
a simultaneous fit of both data sets, the average 
INTEGRAL and average CGROj OSSE spectrum. 
The data are well represented by a cutoff power 
law with individual constant factors in order to 
account for flux variability (xt ~ 1*18). Note that 
the energy range above ~ 100 keV is dominated 
by the OSSE data, as their error bars are signif- 
icantly smaller than the INTEGRAL ones. This 
is because the INTEGRAL exposure time is much 
smaller. The fit shows that the 20-200 keV flux 
varied by a factor of 2 over the individual INTE- 
GRAL and CGRO observations. It is remarkable 
that the same simple model with T = 1.5 and 
E c = 100 keV represents well the data from the 
different CGROj OSSE and INTEGRAL observa- 
tions. 

When applying the PEXRAV model to the 
CGRO data it turns out that a reflection com- 
ponent is not required. This is not surprising, as 
there are no OSSE data below 50 keV where a 
reflection component would have a significant in- 
fluence. The PEXRAV fit to the combined data 
shows that a reflection component is also not re- 
quired by the data. The reduced x 2 is similar to 
that of the cutoff power law, and the fit result is 
consistent with no reflection (R — 0) at a 2 a level. 

This is mainly consistent with previous Bep- 
poSAX observations. Three observations by Bep- 
poSAX , two in 1996 and one in 1999, show that a 
reflection component is not required in most cases 
with R 1996 H = 0 . 01 ,i?i 999 = 0.2, and R 1996 L = 
0.45 (Petrucci et al. 2001). In the latter case the 
fit resulted in a large optical depth r, indicating 
that a reflection component, if present, would be 
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smeared out over a wide energy range making its 
existence more difficult to prove. 

Another study of the high energy spectrum of 
NGC 4151 combining CGRO/OSSE and ASCA 
data has been presented by Zdziarski et al. (2002). 
This work also came to the conclusion that the 
data can be represented by a simple cut-off power 
law with a narrow iron Ka line, adding ther- 
mal bremsstrahlung at low energies, which are 
not covered by our data. This fit resulted in a 
xl = 1.6 and the authors added a broad Gaus- 
sian line in order to account for residuals in the 
5 — 6.5 keV energy band. Zdziarski et al. argue 
that a model with complex absorption and Comp- 
ton reflection component is more likely to give a 
physical meaningful explanation for the observed 
spectrum. They also argue that the absorber in 
NGC 4151 shows evidence for the presence of sev- 
eral absorbing components in the line of sight, 
which cannot be confirmed by our data, mainly 
because of the lack of data below 2 keV. 

Matt et al. (2003) argue that absorption in 
Seyfert 2 galaxies can be caused by several circum- 
nuclear regions. They favour a model with a pos- 
sible temporary switching-off of the nuclear radi- 
ation to explain observed variability. Their model 
predicts an iron Ka equivalent width of — 100 eV 
for a Compton thick case with N h ~ 10 23 cm" 2 , 
consistent with our measurements for NGC 4151. 
This could mean that the iron line could be pro- 
duced in material in the line of sight, with large 
distance to the central engine of the AGN. 

Schurch & Warwick (2002) showed that the 
same ’template’ spectrum matches observations of 
NGC 4151 done by ASCA, BeppoSAX and even 
XMM-Newton (Schurch et al. 2003). Their model 
spectrum includes a strong reflection component 
(R ~ 2 from PEXRAV) as well as an absorption 
component with variable iron abundance. When 
we apply this to the INTEGRAL data with the 
simple cutoff powerlaw and leave the iron abun- 
dance variable, we also detect an enhanced iron 
abundance by a factor of 2.3^’ jp and get. an im- 
proved fit result (xl = 1-33 instead of xl = 1-46). 
In this manner we are able to model the high- 
energy spectrum in a simple way without the need 
of any reflection component. The large error on 
the iron abundance indicates though that the data 
cannot put significant constrains on the true abun- 
dance in the absorbing material. The effect on 


the xl shows that additional components, like a 
variable abundance or a reflection component, can 
improve the fit statistics, but may not lead to im- 
proved insight on the physical conditions in the 
source. 

Most Seyfert galaxies show a softening of the 
X-ray continuum as sources brighten (Markowitz, 
Edelson, & Vaughan 2003; Nandra et al. 1997). 
We did not see significant variations of the spectral 
parameters over the duration of our observations. 
Especially, we did not detect any correlation of 
flux with spectral slope. This can be understood 
when studying the correlation of spectral hard- 
ness ratio with flux as derived from simultaneous 
CGRO /BATSE and RXTE/ ASM data as shown 
in Fig. 6. Here we have plotted the hardness ra- 
tio (20 - 100 keV)/(2 - 12keV) versus the RXTE 
(2 — 12keV) flux based on daily averages from 
1996-2000. An overall correlation is discernible, 
but the scatter is large and thus might not be de- 
tectable over a time span of 400 ksec (as in our 
case), and can even be missed when comparing 
the long-term observations by CGRO/OSSE. 

5. Conclusion 

The results from the INTEGRAL JEM-X2, IS- 
GRI, and SPI data show that the high-energy 
spectrum of NGC 4151 can be approximated by a 
rather simple model. The spectrum above 2 keV is 
described by an absorbed (Affi = 5.1 x 10 22 cm~ 2 ) 
cut-off (Ec = llOkeV) power law (T = 1.5) with a 
narrow iron Ka line (EW = 119 eV). A Compton 
reflection model (PEXRAV) gives a significantly 
better fit result with T = 1.86^’^, a hydrogen 
column density of Affi = 7.0*°‘| x 10 22 cm -2 , a 
folding energy of Ef 0 id = 410*™° keV and a rel- 
ative reflection of R = 1.04*°' 3 [. However, this 
latter result is questionable as the inferred Ej 0 id 
parameter is clearly not constrained by our data. 

The comparison to previous CGRO/OSSE mea- 
surements shows that the simpler cut-off power 
law model is a valid representation of the high 
energy spectrum over decadal timescales, when 
applying a different normalisation to the various 
observations. A significant variation of the spec- 
tral parameters during the INTEGRAL observa- 
tion cannot be detected. The effect of spectral 
softening with higher flux values can only be seen 
in observations covering longer time scales. This is 
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shown in the comparision of CGRO / BATSE (20~ 
200 keV) to RXTE / ASM (2-10 keV) fluxes. 

The INTEGRAL data cannot confirm the ex- 
istence of a high iron overabundance in the ab- 
sorber. 

This is not to say that we claim that Compton 
reflection is not taking place in NGC 4151, but at 
the time of our observations, with the given data 
in the high energy range, we are not able to reach 
a conclusive result. The thermal Comptonization 
is clearly seen in the cutoff of the spectrum, rul- 
ing out non-thermal emission being dominant in 
this source. This is also confirmed by the non- 
detection of emission in the annihilation line in 
NGC 4151. 

The study of NGC 4151 indeed shows that the 
INTEGRAL mission can give valuable insights on 
the characteristics of the high-energy cut-off in 
bright AGN, as has been already shown for the 
Seyfert 2 NGC 4388 (Beckmann et al. 2004a, 
2004b). This shows that AGN can be well de- 
scribed by an optical thick disk with a hot corona. 
We will present a general overview on the prop- 
erties of AGN observed by INTEGRAL in an up- 
coming paper. 

This research has made use of the NASA/IPAC 
Extragalactic Database (NED) which is operated 
by the Jet Propulsion Laboratory and of data ob- 
tained from the High Energy Astrophysics Science 
Archive Research Center (HEASARC), provided 
by NASA’s Goddard Space Flight Center. 
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Fig. 1.— INTEGRAL and RXTE / ASM light 
curves in Crab units. The squares and circles rep- 
resent the 20-100keV ISGRI and 5—20 keV JEM- 
X2 lightcurve, respectively. Fluxes were extracted 
from model fits to the data. The triangles are 
RXTE / ASM 1.5 — 12 keV daily averages. 

NGC 415'i INTEGRAL spectrum 


absorbed cut-off power low with Gaussian tine 



b 


10 'OO 

Channel energy (keV) 

Fig. 2. — Summed spectrum of the INTEGRAL 
JEM-X2, ISGRI, and SPI data. The data are 
well represented by an absorbed exponential cut- 
off power law with Ec — 109 keV and photon in- 
dex T = 1.52. In addition a gaussian iron Kc* line 
with equivalent width EW = 119eV appears in 
the spectrum. 
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Table 1: INTEGRAL observations 


obs. date 
start 

INTEGRAL 

revolution 

exposure 
time [ksec] 

number of 
pointings 

observation 

mode 

22/05/03 

73 

7.3 

3 

5x5 

23/05/03 

74 

170.8 

10 

5x5, staring 

25/05/03 

75 

174.5 

5 

staring 

28/05/03 

76 

55.1 

4 

staring 


Table 2: Spectral fits to average INTEGRAL and OSSE spectra 

Model 

Mission 

Photon 

Eq or kT 

r or R 

Xt ( dof ) 



Index 

[keV] 



Power law 

INTEGRAL 

1 75 + u - u:i 
A *' U -0.02 
9 cc +0.03 

i S ;« 

1 ^ -0.04 



1.90 (207) 

Power law 

CGRO/O SSE 



1.73 (279) 

PL-exp 

INTEGRAL 

109 t\l 


1.46 (189) 

PL-exp 

PL-exp 

CGRO/ OSSE 
INTEGRAL+CGRO 

i 47 + 015 
-0.16 
i 49 +0 04 

A -0.04 

90 +16 
-12 

93 ± 5 5 


0.98 (278) 
1.18 (469) 

Thermal Comptonization a 
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Fig. 3.— INTEGRAL ISGRI (circles, 2000 sec 
binning) and JEM-X2 (triangles, 1000 sec binning) 
lightcurve in the same energy bands as in Fig. 1 . 
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Fig. 4. — INTEGRAL/ OMC light curve. Error 
bars are lcr values. 



KiGC 4151 combined INTEGRAL and OSSE spectrum 
PEXRAV model; photon index 1.6, fold* 123 keV, refl. 0,41 



Fig. 5.— INTEGRAL and CGRO / OSSE simul- 
taneous data fit. The OSSE data represent the 
average of the NGC 4151 spectra. All data are 
well represented by a PEXRAV model, with the 
individual OSSE data scaling with 0.6 ... 1.2 ver- 
sus the INTEGRAL/ JEM-X2 data. 
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Fig. 6. — Spectral hardness derived from the 
BATSE and RXTE data, versus the RXTE / ASM 
flux. A softer spectrum is correlated with a higher 
flux. The error bars on the lower left show the 
typical uncertainties of the measurements. 
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